Parametric investigation of nonlinear fluctuations in a dc glow discharge plasma 
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Glow discharge plasmas exhibit various types of self excited oscillations for different initial condi- 
tions like discharge voltages and filling pressures. The behavior of such oscillations associated with 
the anode glow have been investigated using nonlinear techniques like correlation dimension, largest 
Lyapunov exponent etc. It is seen that these oscillations go to an ordered state from a chaotic state 
with increase in input energy i.e. with discharge voltages implying occurrence of inverse bifurca- 
tions. These results are different from the other observations wherein the fluctuations have been 
observed to go from ordered to chaotic state. 

PACS numbers: 



Chaos and different route to chaos have been 
studied in many glow discharge plasma system. 
In the earlier experiments which are mainly per- 
formed in the parallel plate electrode systems, 
it is observed that the system becomes chaotic 
with increasing input energy. But in cylindri- 
cal system, we have observed an opposite phe- 
nomena i.e. chaotic to stable state transitions 
with increase in input energy, where the energy 
is given the form of the discharge voltages. This 
system has extensive use in the dusty plasma ex- 
periments. We provide a qualitative and quan- 
titative description, based on ideas and methods 
from nonlinear time series analysis, of transitions 
in these dynamics. We also provide a plausible 
model for this phenomenon. We hope this work 
can help improve the fundamental understanding 
of such inverse phenomenon. 



I. INTRODUCTION 

Plasma is a typical complex medium exhibiting a wide 
variety of nonlinear phenomena such as self oscillations, 
chaos, intermittency etc [H, [H H, H[ . These systems ex- 
hibit chaotic or regular behavior depending on the con- 
trolling parameters like discharge current, voltage, pres- 
sure etc, and the transition from one state to another 
takesplace for a change in any one of the parame- 
ters [J H, H, 0, [1] • In this paper we report on the ob- 
servations of the chaotic and the regular behaviors in the 
floating potential fluctuations at different discharge volt- 
ages in a concentric cylindrical electrode system with a 
hollow cathode. In almost all the earlier experiments, the 
system goes from order state to chaotic state, whereas in 
ours the reverse is seen, thus indicating signature of in- 
verse bifurcation for the first time in a glow discharge 



plasma. The fluctuations are observed when a bright 
glow is formed around the anode. Nonlinear fluctuation 
analysis has been carried by estimating the correlation 
dimension and Lyapunov exponent in addition to power 
spectral analysis. The correlation dimension and Lya- 
punov exponents were estimated using the Grassberger- 
Procaccia algorithm 0, the Rosenstein's technique algo- 
rithm [I(J. The latter was chosen for the Lyapunov ex- 
ponent since it is efficient for small data lengths [ljl [H| • 



II. EXPERIMENTAL SETUP 

The experiments were performed in a hollow cathode 
dc glow discharge plasma. The schematic diagram is 
shown in Fig Q] A hollow S.S. tube of diameter (cf>) w 
45 mm was the cathode and a central rod of <f> ~ 1.6 mm 
was the anode. The whole assembly was mounted inside 
a vacuum chamber and was pumped down to about 0.001 
mbar using a rotary pump. The chamber was filled with 
argon gas up to a desired value of neutral pressure, by a 
needle valve, and then a discharge was struck by a dc volt- 
age (which we called as discharge voltage or DV), which 
could be varied in the range of 0—1000 V. The electro- 
static floating potential fluctuations were measured using 
a Langmuir probe of <f> — 0.5 mm and 2 mm long, approx- 
imately at the mid position of the anode and the cathode. 
A black and white CCD camera has been used to measure 
the size of the anode glow. The measured plasma density 
and the electron temperature were about 10 7 cm~ 3 and 
3—4 eV respectively. The corresponding electron plasma 
frequency (f pe ), and the ion plasma frequency (f P i) were 
about 28 MHz, and 105 kHz respectively. 



III. RESULTS 



A. Model for anode glow 
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By varying the neutral pressure, we observed the 
breakdown of the gas at different voltages as shown in 
Fig [21 The breakdown voltage (V& r ) initially decreases 
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FIG. 1: Schematic diagram of the cylindrical electrode sys- 
tem of the glow discharge plasma. The probe is placed at a 
distance / ~ 12.5 mm from the anode. 



with increase in pd, (where p and d are the filling pres- 
sure and radius of the cathode respectively) and then 
begins to increase with pd after going through a mini- 
mum value resembling a typical paschen curve [l2j]. But 
the characteristics of the discharges are very different at 
low and high pressures. At higher pressures especially for 
pd > 20 mbar-mm an anode glow was observed as shown 
in fig 3. 
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FIG. 2: Vbr vs pd (Paschen curve) for our experimental sys- 
tem. The minimum occurs at (1.69 mbar-mm, 251 V). 



Figs [3£a) shows that the glow with largest size, ap- 
pears when the discharge is struck at a typical pressure 
of 0.95 mbar and its size decreases with increase in the 
DV until it finally disappears [Figs [3Ja)-E[h)]. From 
the CCD image analysis the annular radius of the glow 
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FIG. 3: Evolution of glow size of the anode glow with increas- 
ing DV (a) - (h). 



around the anode, was estimated to be w 1.3 mm at 
the beginning of the discharge [Fig^a)] and reduced to 
w 0.32 mm [Fig [3Tf)1 . We estimated the thickness (5) 
using the relation [13[ 



5 w 3.7 x 10" 



. kT 



(i) 



where k, T, u,-, and P are Boltzmann constant, room 
temperature in Kelvin scale, ionization cross section and 
presure in mbar respectively. Ci depends upon the energy 
of the electron and at 15.76 eV, a t is « 2 x 10~ 18 cm 2 [13 ]. 
The estimated thickness of the anode glow for P = 0.95 
mbar; T = 300 K, is « 0.81 mm which is within the 
range of the thickness estimated from the image shown 
in Fig [31 

The anode glow appears when the discharge current 
is too low to sustain the discharge [HI, [l5|. This glow 
supplies energy, in the form of positive potential gradient, 
to the electrons so that additional ionization takes place 
and hence discharge current increases [H, Ell- In the 
present experiment, at the initial stages of discharges the 
electrons do not have sufficient energy near the anode to 
reach it, and hence the anode glow occurs. But, when 
the electrons have sufficient energy to reach the anode 
the glow disappears. 
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FIG. 4: Spacial variation of the plasma floating potential: (a) 
without anode glow; (b) with anode glow. Vertical dotted 
and dashed lines are the cathode positions (C) and anode 
position (A) respectively. Horizontal dash dot line is the zero 

level. in the curve near the anode and the cathode are 

the extrapolation of the potential. 



FiglUa) andHfb) show the tentative model of floating 
potential profile for the present experimental arrange- 
ment. It shows that the floating potential increases 
sharply in the sheath region near the cathode and then 
decreases slowly up to the anode. As the measurement 
near the anode and the cathode was not possible, the 
extrapolated potential profile is shown by dash lines. In 
the the presence of the anode glow, the potential near 
the anode get modified. Experimentally it is observed 
that the anode glow is separated from rest of the plasma 
by a double layer [l3|, [3, U3] and the height of the 
potential of the core of this glow is of the order of the ion- 
ization potential [l6j], which for argon is about 15.76 eV. 
As the anode is grounded, the modified potential profile 
in this case, over the profile shown in Fig0Ja), will be as 
shown in Fig0|b). As the width of the anode glow de- 
creases with increase in the DVs [Fig [3] the hump shown 
in Fig HJb) will also shrink accordingly. This shrinkage 
happens because as the electron energy increases with the 
DVs, the necessity of the additional ionization to main- 
tain the discharge decreases. 



B. Analysis of the floating potential fluctuations 

An interesting feature associated with the anode glow 
was the different types of oscillations in the floating po- 
tential at different pressures. 

The behavior of the floating potential fluctuations for 
pressures less than the Paschen minimum, was different 
from that at higher pressures. At low pressures (LHS of 
Paschen minimum) the fluctuations exhibited Self Orga- 
nized Criticality (SOC) like behavior but at higher 
pressures (RHS of the Paschen minimum), we observed 
fluctuations in the floating potential with the formation 
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FIG. 5: Sequential change in the raw signal (normalized) at 
0.89 mbar for different voltages: (a) 288 V; (b) 291 V; (c) 295 
V; (d) 301 V; (e) 304 V; (f) 307 V; (g) 327 V; (h) 385 V; (i) 
466 V. All y-axes range form -1 to 1. 



of the anode glow (fire ball, similar to Ref [13|, Hsj ) 
which also disappears with the disappearance of the glow 
[Fig . 

We carried out a detailed analysis of the fluctuations 
for three typical pressures which are presented here. At 
about 0.89 mbar (pd ~ 20.02 mbar-mm), the discharge 
was initiated at w 288 V and an anode glow was observed 
similar to Fig [3Ja) . Simultaneously irregular relaxation 
type of oscillations in the floating potential were observed 
[Fig Eta)]- Increasing the DVs, led to an increase both 
in the amplitude and the time period of the oscillations 
[Figs O (b)-(Sfi)]. The regularity of the oscillation also 
increases with the DVs. Around 509 V both the anode 
glow and the fluctuations disappeared simultaneously. At 
0.95 mbar (pd w 21.37 mbar-mm) the fluctuations are 
observed to be more random [Fig [TJJa) — (i)] than at 0.89 
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FIG. 6: Sequential change in the raw signal (normalized) at 
0.95 mbar for different voltages: (a) 283 V; (b) 284 V; (c) 286 
V; (d) 288 V; (e) 289 V; (f) 290 V; (g) 291 V; (h) 292 V; (i) 
293 V. All y-axes range form -1 to 1. 

mbar. With increasing DV, the final form of the fluctua- 
tions before their disappearance was the relaxation type 
of oscillation as shown in Figs EJh) — Efi) . Increasing the 
pressure causes more randomness in the signals as seen 
in Fig[7{a) — (i) at 1.0 mbar (pd 22.5 mbar-mm). We 
observed the relaxation type of oscillations at DV « 313 
V just before the glow an the fluctuation disappear. 



Figs EKl}, {II} and {III} show the power spectrum 
calculated for the signals at three different experimen- 
tal pressures shown in Figs [5J H] and [7] respectively. A 
broad range of frequencies (< ion plasma frequency) is 
observed at the initial stage of the discharge, at pres- 
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FIG. 7: Sequential change in the raw signal (normalized) at 
lmbar for different voltages: (a) 293 V; (b) 296 V; (c) 298 V; 
(d) 299 V; (e) 300 V; (f) 305 V; (g) 308 V; (h) 310 V; (i) 312 
V. All y-axes range form -1 to 1. 



sures 0.89 mbar [FigEp} (a)], 0.95 mbar [FigHJII} (a)] 
and 1.0 mbar [Fig[8]{III} (a) — (b)]. The figures also show 
that with increasing DV, number of peaks present in the 
power spectrum, decreases. At higher DV i.e. at the end 
of the sequence of the power spectra we observed one ma- 
jor peak with its harmonic [Figs [8]{I}(f) — (i) and III(i)] 
or two peaks [Figs [5][II}(i)] just before the fluctuations 
cease. The power spectra clearly show that the system 
goes from an irregular to a more or less regular state with 
increasing DV. Another interesting feature that has been 
observed for all the three pressures is that the dominant 
frequency i.e. the frequency which has maximum power 
decreases with increase in the DVs [Figs ETa) — (c)]. On 
the contrary, in many experiments [l9j, [20j, l2l| , the domi- 
nant frequency is observed to increase with the DV. The 
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FIG. 8: Power spectrum of floating potential fluctuations at 
initial discharge voltages at filling pressures: (a) 0.89 mbar, 
(b) 0.95 mbar and (c)l mbar. 
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FIG. 9: Dominant frequency vs DV at the three experimental 
pressure: (a) 0.89 mbar; (b) 0.95 mbar; (c) 1.0 mbar. 



basic difference between the present experiment and the 
experiments performed in Ref [lfj, |2Cj, |2l| was in the elec- 
trode configuration and biasing arrangement. In those 
experiments [T^ H(| [U the system was planar type is 
which the anode voltage was increased keeping the cath- 
ode at ground whereas ours is a cylindrical system in 
which we increased the negative voltage on the cathode 
keeping anode at ground. In those systems electrons orig- 
inating at cathode get lost at the anode. In our case the 
electrons can undergo multiple reflections and hence in- 
creases the ionization length. The fluctuations are prob- 
ably due to ion acoustic instability driven by electron 
beam plasma interaction and also from electron and ion 
trapping in the potential wells, associated with the anode 
glow [Fig [4]. 

Ions produced inside the anode glow due to collisions 
of the accelerated electrons across the hump with the 
neutrals, makes the glow unstable [l5|, which is proba- 
bly responsible for the relaxation or the random oscil- 
lations [l^, H(|. An estimate of the frequency of these 
instabilities can be obtained from the ion transit time in 



electrode distance. The estimated ion transit frequency 
(-) for our experimental system is ~ 19 kHz which agrees 
well with the frequency of the relaxation oscillations of 
the floating potential shown in Figs 0[I}(f) — (i),{II}(i) 
and {III}(i). The higher frequency oscillations could be 
due to ion acoustic instabilities since the conditions are 
quite conducive to excite these low frequency instabil- 
ities. In the next section we have presented nonlinear 
analysis of the fluctuations. 

The presence of relaxation oscillations have been at- 
tributed to the formation of highly nonlinear structures 
like double layers [12]. We therefore estimated the cor- 
relation dimension (D corr ) and the -l-ve Lyapunov ex- 
ponent (Al) of all the signals. In this experiment, it is 
observed that the nature of the fluctuations of the poten- 
tial did not vary for almost the whole day as long as the 
controlling parameters are kept constant, so the analyzed 
signals can conveniently be taken to be stationary. 

The correlation dimension (D corr ) has been calculated 
using well-known Grassberger-Procaccia techniques [^, 
[23j . where correlation sum (C(r,m)) scales with scale 
length (r) as C(r, to) = r^ orr . A typical plot of In C(r, to) 
vs In r has been shown in Fig [TO] for the embedding di- 
mension (m), in the range of 2 to 10, from which we have 
estimated the correlation dimension (D corr ) for a typical 
signal (FigE^a)) at DV w 283 V, at 0.95 mbar. The scal- 
ing region had to be chosen carefully, since for too small 
scale lengths the correlation sum is heavily distorted by 
noise and the higher scale lengths are limited by attrac- 
tor dimensions. From the above plot the correlation sum 
exhibits a power law behavior within a certain range of 
r as shown by the vertical dotted line. In C(r, m) vs lnr 
plots are almost parallel for higher m (i.e. m = 7 — 10) 

and the corresponding best fit has been shown by 

line. D corr vs m is also shown in the inset of Fig [TO] It 
shows that D corr becomes constant at higher m and this 
constant value of D corr is the correlation dimension of 
that particular signal and for all of our data D corr > 3.8 
to begin with. The D corr at 0.89, 0.95 and 1 mbar 
for different DVs have been shown by open circle(o) in 
FigsHlTal-fTlTc) respectively. It is observed that there is 
a decreasing tendency of D corr except for some interme- 
diate values of DV at higher pressures, where it is seen to 
increase and then decrease again. Since D corr is a mea- 
sure of the complexity of the system it is likely that the 
system complexity increases at those intermediate DVs. 
In our experiment initially the system is in a complex 
state as D corr for all the three pressures is high and de- 
creases with increase in DVs [FigfTTTa) — (c)]. D corr ~ 1 
just before the system reaches stable state finally, is the 
indicator of the periodic state and these periodic nature 
are also prominent from the raw data and the power spec- 
trum [Figs [7]and Fig|S][I}-{III}]. D corr shows that 
the system stabilizes itself with increase in DVs. 

The presence of a -l-ve Lyapunov exponent (Xl) is the 
most reliable signature of the chaotic dynamics and it 
is estimated using the Rosestein algorithm [l0( which is 
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FIG. 10: Effect of embedding dimension on correlation sums. 
The scaling region is shown by two vertical lines. The best 
fitting is shown by line. 
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FIG. 11: D corr vs discharge voltages. +ve Xl has been shown 
by +ve sign, at (a) 0.89 mbar; (b) 0.95 mbar; (c) 1 mbar. 



especially suited for small data lengths. Xl can be es- 
timated from the slope of the plot < \ndj(t) > vs time 
(iAt) as shown in Fig [121 The figure shows < lndj(t) > 
vs iAt for m=2 to 10. A clear scaling region is seen at 
a higher m shown by line. The positive Xl has been 



TABLE I: Error in the estimation of the correlation dimension 
and the +ve Lyapunov exponent 
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FIG. 12: Average In dj (t) for different embedding dimensions. 
The best region for Xl has been shown by dotted line. 

identified at 0.95 and 1 mbar for some DVs and they are 
shown in Fig HHb) and [TTT c) by +ve sign respectively. 
Figs fTTT b) and (c) show that Xl becomes positive for 283, 
284, and 290 V at 0.95 mbar and for 293, 296, 300 and 
305 V at 1 mbar respectively. Though D corr quantifies 
the complexity present in the system, it does not guar- 
antee the presence of chaos, which is determined by Xl ■ 
At 0.89 mbar initially we have a high D corr [FigfTTTa)]. 
but the Xl is not positive in this case, implying that the 
system is not in chaotic state. At higher pressures we 
find that in general Xl is +ve D corr > 3, suggesting a 
low dimensional chaos. 

D corr and A^ have been estimated for different data 
lengths to check the error on the estimations of the D corr 
and Xl- The estimated D corr and Xl at different data 
lengths have been shown in Table |TJ Both and A l and 
Dcorr tend to show stable results for higher data lengths 
(larger than 1000) as seen from Table HI 

The evolution of the floating potential fluctuations 
have also been shown using the bifurcation diagram in 
Fig [HI At 0.89 mbar the system reaches a stable state 
through period subtraction [Fig [lljTa')] with increasing 
DV. No chaos was observed in this case. For 0.95 mbar 
a chaotic region between DVs 283—284 V and an in- 
termediate chaotic state have been observed at 290 V 
[Fig [T37 b)]. We have two chaotic regions, one between 
288-296 V and the other between 300-305 V at pres- 
sure 1 mbar [Fig [PBT c')]. The figures show that inverse 
bifurcation probably takes place in the system and this 
has been observed for the first time in a glow discharge 
system. 

In order to validate the above observations, we carried 
out the surrogate data analysis to detect the stochastic- 
ity or determinism in the system 0, Hfl HE| • We have 
employed the correlation dimension as the test statis- 
tics for this method. The surrogate data has been gen- 
erated by Phase Shuffled method, in which the phases 
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FIG. 13: (a), (b), and (c) are the bifurcation diagram with 
DVs at pressures (a) 0.89, (b) 0.95, and (c) 1 mbar respec- 
tively. 
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FIG. 14: Surrogate data using phase randomizing method: 
(a) original data; (b) surrogate data; (c) Correlation dimen- 
sions vs m for original (A) and its surrogate data (©). 



of the signals are randomized by shuffling the fourier 



phases [2J, [25|, l27| , and hence the power spectrum (lin- 
ear structure) is p reserved, but the nonlinear structures 
are destroyed [25| . As seen in Fig H4T c) The D corr for 
the original data saturates at higher m, whereas for the 
surrogate data it increases with m as expected, since the 
surrogate data is supposed to be random in nature, and 
its D corr should be infinite |27| . The saturation of D corr 
in our case indicates that the system dynamics is deter- 
ministic in nature. 

IV. CONCLUSIONS 

Glow discharges are simple systems but exhibit exotic 
features depending on the configuration, initial parame- 
ters etc. Since various configurations are used in different 
applications like dusty plasma, plasma processing etc. it 
is necessary to understand the plasma dynamics of the 
these systems as much as possible. However the under- 
standing of the complexities in the plasma dynamics is 
quite a challenging job as they arise from many degrees 
of freedom like different sources of free energy, different 
types of wave particle interaction and many other insta- 
bilities. In our present work nonlinear time series analysis 
has been used to quantify and differentiate complex and 
coherent processes at different parametric conditions. We 
also observed that the chaotic state driven by a driving 
force (DV) tends to be stabilized by some type of inverse 
bifurcation. This is also evident from the estimated D corr 
and Xl at different DVs. It is not fully understood as to 
why relaxation oscillations are being observed only at 
0.89 mbar and not at 0.95 and 1.0 mar though there is 
not much difference in the three pressures. One probably 
has to carry out a detailed simulation to investigate the 
evolution of the fluctuations at different discharge volt- 
ages. 
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